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1 Introduction 

In a recent paper [Ij (hereafter referred to as ACV) a simplified non-perturbative treat- 
ment of transplanckian scattering has been proposed. It is based on an approximate re- 
summation of the semi-classical corrections to the leading eikonal approximation. These 
were identified long ago [2] as coming from tree diagrams where the fast initial particles 
act as external sources for gravitons. 

Resumming tree diagrams amounts to solving the equations of motion of an effective 
action proposed quite sometime ago [H H] . Taking the high-energy limit simplifies consid- 
erably the longitudinal dynamics which, in a certain approximation (neglect of so-called 
"rescattering diagrams"), can be factored out leaving behind an effective dynamics in the 
"transverse" [D — 2)-dimensional space. In pp (and here) the case oi D = 4, hence of an 
effective two-dimensional theory, was (will be) considered. 

In [1] one further simplification consisted of neglecting a physical polarization of the 
radiated gravitons, the one suffering from infrared problems. Finally, some azimuthal 
averaging was used in order to be able to treat the problem both analytically and nu- 
merically. In spite of these approximations, the resulting S-matrix looked able to capture 
the essential physics of the gravitational collapse problem. In particular, it showed the 
existence of a critical impact parameter be below which a new absorption (and thus the 
opening of some new channels) occurs. That critical value, possibly signalling a transition 
between what we may call a dispersive (D) and a black-hole (BH) phase, turned out to be 
in good agreement with the one based on closed-trapped-surface (CTS) criteria [5], [S]- 
A more recent investigation [7] managed to solve numerically the full partial differential 
equations (PDE) (i.e. without using the azimuthal average approximation) and largely 
confirmed the qualitative outcome of the simplified ACV approach. 

In this paper we consider more systematically the case of axisymmetric extended 
colliding sources/beams. The motivations for considering this particular case are several: 
on the one hand the PDE's reduce to ODE's making the problem an affordable one by 
analytic (or by much easier numerical) techniques, without having to make the azimuthal 
averaging approximation. At the same time, as we will show, the IR-sensitive polarization 
is simply not produced in the axisymmetic case. Last but not least, we are able to 
introduce a vaste number of initial states, by playing at will with the many (shape and 
intensity) parameters chracterizing the sources and check for the existence of critical 
surfaces in this multidimensional space. The results can then be compared with those 
coming from CTS criteria [6] and, hopefully in the near future, will be tested against 
analytic [8J as well as numerical GR calculations [9]. 

The paper is organized as follows: in section 2 we make some general considerations 
on the axisymmetric case and prove a one-way relation between the CTS criterion of [6J 
and the criticality condition in the ACV system. In section 3 we give some examples of 
interesting extended sources while in section 4 we determine numerically the critical lines 
in the sources' parameter space. We also compare those lines with the bounds coming 
from the (much simpler) CTS conditions and find very good agreement. In section 5 we 
study the behaviour of the elastic S-matrix in the vicinity of the critical points and point 



out a possible connection with Choptuik's scaling ^U\ near critical collapse. Section 6 
briefly summarizes our results and future prospects. 

2 The axisymmetric case: general considerations 

Our starting point is the effective two-dimensional action of [Ij (see their equation (2.22)): 
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-V'^n = V'^aV^a-ViVjaViVja, (1) 

where a, a and are three real fields representing the two longitudinal and the (IR-safe) 
transverse component of the gravitational field, respectively. Equation ([1]) can be easily 
generalized in order to deal with two extended sources: 
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where the center of mass energy ^/s provides the overall normalization factor 2txGs = 
^i?^, while the two sources s{x),s{x) are normalized by J d'^x s{x) = J d'^x s{x) = 1. 

Let us now specialize to the case of two extended axisymmetric sources moving in 
opposite direction with the speed of light and undergoing a central collision. Using the 
conventions of [1] we will denote by Eiiri) (in the following i = 1,2 will represent unbarred 
and barred fields/sources respectively) the energy carried by the ith beam below r = ri 
and define Ri{r) = AGEiir). Let us also assume that the two sources have finite support 
so that Ri{r) = Ri{oo) = Ri for r > L^. By going to the overall center of mass, we may 
always choose Ri = R = 2Gy/s. 

2.1 Simplifications 

One advantage of considering the axisymmetric case is that there is simply no dependence 
of the physics upon the azimuthal angle, hence no need to take averages over it. This is 
a useful technical simplification that allows us to reduce the problem to solving ODE. 

The second more important advantage comes from the observation that the IR-singular 
"LT" graviton polarization is not produced in that case. Thus the problem is completely 
IR-finite even in D = 4. In order to see this, let us recall from [Ij that the LT polarization 
is produced with an amplitude proportional to sin 6*12 cos 6*12. In the notations of yj: 

Alt = A^^e^j^j. ~ k~^ sin 6112 cos 6112 , (3) 

where ^12 is the angle between the two transverse momenta fei, fc2 that combine to give a 
physical graviton of momentum k. The angular factor can be expressed in terms of the 



momenta as: 

klkl sin 6112 cos 6112 = eijk\k\klk2 ■ (4) 

Going over to position space this becomes a differential operator acting on a and a: 

eijV'V^aS/^V^a . (5) 

Given that in the axisymmetric case a and a depend only on x'^ it is easy to see that 
the above expression vanishes as a result of the antisymmetry of eij. 

2.2 Axisymmetric action and equations of motion 

It is straightforward to rewrite the action ([2]) for the axisymmetric case as a one dimen- 
sional integral over the variable r^ = x^ = t. Using J (Px = tc j^ dt we find: 
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dt{l-p)\ (6) 



(27ri?)2 
where a dot means d/dt and, as in [Ij, we have introduced the field: 

p = t(l- (27ri?)20) . (7) 



Integrating by parts and using tt J^ dt'si{t') = Ri(t)/R we arrive at the following conve- 
nient form of the action: 
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The equations of motion that follow from ([8]) read: 
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and therefore reduce to a closed 2nd order equation for p. We want to look for solutions 
of that equation with the following boundary conditions [Ij: 

p(0) = , p{t)^t as t^oo. (10) 

Given the finite support of the sources the latter condition can be replaced by the re- 
quirement: 

p = V^l - R^/p for y/t> Max{Li,L2).. (11) 



For given source profiles Ri{t) a possible strategy for solving the problem is to reduce 
it to a first order system: 
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with initial conditions 

p(0) = , a(0) = ao, , (13) 

and to find a cto such that cr(Maa;(Li, L2)) = 1. For sufficiently small (large) Ri/Li 
one expects that the latter condition can (cannot) be imposed on real-valued solutions. 
Thus, in general, there should be some critical values for those parameters separating two 
different regimes. 

2.3 CTS-criteria and critical points: a general result 

In the general axisymmetric case, one can construct explicitly a minimal CTS [6j provided 
that an r^ exists such that (see eq. (4.4) of [6J for D = A): 

Riir,)R2ir,) = rl (14) 

We will now argue that such a condition implies the absence of real solutions to eqns. 
([9]) with p(0) = 0. Proof: Let us first note that, because of flT2l) and the fact that the Ri 
are non-decreasing functions of r, the quantity a, as well as p, are increasing functions of 
t. Therefore, for any t: 



a(t)<a(oo) = l , i.e. pW < W 1 " ^^^^^^ • (15) 

Assuming that the KV criterion (fT4|) can be met let us write: 



p(0) = p{Q - I dt'p{t') > p{Q - t,p{Q > p{Q - U 1 - -^ , (16) 
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where we have used eqs. (1141) and flTSj) . At this point it is easy to check that the rhs of 
flT6|l cannot vanish for any (positive) value of p{tc) thus proving that we cannot impose 
the condition p(0) = when the criterion flT4l) is satisfied. 

3 Examples of source profiles 

3.1 Two identical two- parameter sources 

A class of finite-size sources is characterized by the following profiles: 

d 
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where, without lack of generahty, we have fixed the sizes of the two beams to be 1. One 
can easily verify that these sources satisfy : 

TT / dt's(t') = R(t)/R = -T^r , TT / s(t)dt = 1 . (18) 

Jo {d+il-d)t^f^ Jo 

Inserting these expressions in our differential equation (fT2|) leads to the equations: 
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R\t) = ^ , ;^ ^ ,,,, , t<i, (19) 
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while for t > 1 we should impose a = 1. 

This system can be easily studied numerically. In particular one can find how the 
critical value of i?, Re-, depends on the parameter d. This parameter gives (for fixed total 
extent L = 1) the "shape" of the extended sources: d = 1 corresponds to the case of 
constant-density sources considered in [Tj , (i < 1 to sources roughly concentrated around 
r ~ ydL^ and d> 1 to sources peaked around r = L = 1. 

3.2 Central scattering of a particle off a ring 

In this, very asymmetric case Ri = R while R2 = RQ{t — IP'). The notation anticipates 
that we can view this case as an approximation to the scattering of two particles at impact 
parameter h. And indeed, from ([9]), we recover in this case the approximation used in [I] 
to describe the latter process i.e. 



lRi{t)R2{t) _ IR^ 
2 ^2 - 2^ 



p='^tliy:^>^P^ = ttl_Qit-h^). (20) 



If we require p(0) = this equation leads to the condition on p(6^): 

p(6^) = 6V(fe^)=6%/l-^, (21) 



which has a real solution only if ^ < (y)c = 2^/23 ^^^ ~ 0.62. Such a result has to 
be compared with the upper limit given by (TT4|) which, in this case, simply becomes 
{R/hfJ' < 1. 

It is interesting to notice that exactly the same (y)c will apply to the situation in 
which the point-like source is replaced by an arbitrary source "contained" inside the ring- 
shaped one. Physically this makes sense since, by Gauss' theorem, the compact source 
should propagate undisturbed while the more extended source is only affected by the total 
energy of the more compact one. 



3.3 Gaussian sources 

Point-like sources are difficult to deal with numerically (especially in momentum space). 
Therefore we also introduce Gaussian-smeared versions of the above point and ring-like 
sources. 

s,{x) = i.exp (-^)0 iL,-r), s,{x) = i.exp (-^^^^) ©(^^ - r), 



M - 2.aM - exp {0) ) , M ^ 2. (.^(exp (^) - 1) + -^./fErf-^ 
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When a — > (oo) the problem reduces to the one of the point-ring (two homogeneous 
beams) case. 

Another interesting example is the one of gaussian sources concentrated at r = 0. 
They correspond to: 

1 f t \ RAt) f t \ 

sm = ^-^exp ^-^ J , -^ = l-exp ^-^J . (23) 

4 Numerical solutions and critical lines 

Solving Eqs.(l9]) serves two purposes: a) to determine a critical line which separates the 
dispersive (D) and black-hole (BH) phases, and b) to obtain the actual solution, for a given 
source parameters, and to calculate various on shell observables. Let us illustrate the ffist 
procedure in the case of the scattering of the two identical beams with the sources fITTj) . 
There are two independent parameters: d - specifies the shape of the energy distribution, 
and R/L the ratio of the Schwarzschild radius of a source to its maximal extent L. We 
shall work in units of L. As already outlined in Sect. 2.1 we would like to find an initial 
value (To, in (fT3|) . such that the corresponding solution satisfies cr(l) = 1, matching the 
t > L^ behaviour (TTT]) . We shall refer to such a solution as the maximal solution. 

To this end, we examine the dependence cr(l) on ctq, at given {d, R), cf. Figure [1], and 
identify the region, in {d, R), where the equation a{l, ctq) = 1 has a real solution. It turns 
out that when we increase R, at fixed d, the minimum in FigUl moves upwards giving 
rise to three cases: i) for R < Ri there is only one solution, ii) for Ri < R < R^ there are 
two (real) solutions, and iii) for R^ < R there are no real solutions. We then determine 
the critical values Re for a range of d thus obtaining the critical line Rc{d) = {R/L)c{d) 
in the [d, R/L) space. All numerics is done with the aid of Mathematica. 

The final result, cf. Fig|21 con&ms well our expectations. In the limit of the homoge- 
neous sources {d = 1) we recover R^ = 0.47 - the value already quoted in [IJ. For d — > oo 
one expects to reach the kinematics of the particle-ring scattering, with Re = 0.62, and 
indeed the critical line is consistent with this prediction. In Fig|2] we also plot the upper 
bound on R^ coming from the GTS criterion (fl^ . It is easy to check that this reads: 

rCts ^ (4^(1 _ ^))i/4 0(^/2 -d) + e{d - 1/2) . (24) 



a(l) 1.1 - 




Figure 1: Looking for the maximal solution of Eqs. fll2l) . a{l) 
initial value cr(0) = ctq (in the specific case d = 1, R/L = 0.46) 



1, as a function of the 



Log Re 



Log d 



Figure 2: The critical (solid) line in the {R,d) plane having set L = 1. We also show 
(dashed hne) the upper bound on Re from the CTS criterion flT^ . The BH phase is above 
the solid line. 



We have also checked that, to a good approximation, our Re roughly varies as d^^^ for 
small d in nice (though not necessary) agreement with flM|) . 

As a second example consider the scattering of two Gaussian sources, concentrated at 
a; = with transverse sizes Li, L2, as in fl2^ . Strictly speaking, the sources do not have a 
finite support, however they vanish quickly for t ^ Li. Therefore we replace the condition 
cr(l) = 1 with cr(#Max(Li, L2)) = 1. In practice taking # = 10 is more than sufficient. 
This time we set i? = 1 and examine the transition in the (Li,L2) plane. The result is 
shown in Fig. 3. Obviously, the critical line is symmetric with respect to Li ^^ L2. It's 
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shape corresponds roughly to a straight hne Li + L2 = const (as clearly visible in a linear 
plot not shown here), suggesting that this variable controls the effective concentration of 
the total energy. Once more we can ask how this critical line compares with the lower 
bound one would find using the GTS criterion flT^ . This can be easily done numerically 
and the result is shown again in Fig.3. Amusingly, also the GTS criterion gives a curve 
which, in the central region Li ~ L2, roughly corresponds to Li + L2 = const. Besides 
a typical factor 2 discrepancy (but remember: GTS criteria only give bounds!), the two 
curves start to diverge when the two Li are very different. In fact, while the GTS-curve 
goes to the point Li = 0,^2 = -j^, our result indicates that this point moves to infinity 
as Li -^ 0, i.e. that one is always in the BH regime in that limit. 



Log L2 




Figure 3: Two gaussian sources around the origin. The BH phase lies below the critical 
(solid) line in the (Li, L2) plane. We also show (dashed line) the lower bound on the curve 
from the GTS criterion (fT4l) . 



All the above calculations required solving Eqs.( IT2|) for t < 1. Next we look for the 
maximal solution in the larger domain t < t^ax- To this end it suffices to use the value (Tq 
found above and extend the profile Riit) beyond L^: Riit) = R,t > L"^. The rest is again 
done by Mathematica. As an example, we show in FigH] the derivative of the solution 
p(t). It interpolates smoothly around t = 1 and, as expected, tends to 1 at large t. We 
have also verified that for t > 1, p{t) obtained above is identical to the analytic solution 
for the constant profiles 

Panit) = R^F-' [Fipo/R") + t/R^ - to/R"] , (25) 

with po = p(to), to > 1, and 

F{x) = ^x{x - 1) + log{^/x + ^x - 1). (26) 




Figure 4: The maximal solution, in the dispersive phase, also for t > 1 (two identical 
sources, eq. ([17]), with d = 1, R/L = 0.44). 



5 Behaviour near the critical point 



An interesting quantity to compute is obviously the on-shell action, i.e. the value of the 
action on the e.o.m. as a function of the external parameters (the profiles of the external 
sources). In the multidimensional parameter space we expect to find, as in the collapse 
problem, a critical surface of co-dimension 1 where the on-shell action itself develops 
some branch-point singularity. Typically, on the other side of that critical surface the 
action acquires an additional imaginary part (on top of the one due to graviton emission) 
providing an additional absorption of the elastic amplitude |1] . 

This extra absorption could possibly be interpreted as the opening-up of new channels 
e.g. those corresponding to the formation of black holes behaving as stable particles at 
leading order in h. The absorption itself is expected to be related (up to some numerical 
factor) to the number of new channels that are opening up beyond criticality and to 
provide therefore information on the entropy of the black holes that are supposedly formed. 
Under these assumptions a connection can possibly be established between the critical 
exponents of the action and those of Choptuik's scaling ^0\. 

On the equations of motion the action ([8]) immediately simplifies to read: 



A' 
h 



4/2, 



dt 



[l - pf - -R,{t)R2{t) 
P 



(27) 



while using the second of eq ([9]) the above expression can also be written in the form: 



A eon 



1 



dt [(1 - pf - 2p-p\ 



.^Jdt[3{l-pf-2tp] 



1 



dt 



3{l-pf^2tp-2-Ut-p){l-p)] 



dt' 



(28) 
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Here we have used the fact that the last surface term vanishes because of the properties 
of the solution at t = and t = oo. 

We would like to study the action near the critical point where we expect to find 
a branch point behaviour. Unfortunately, the action itself is IR-singular: the infinite 
Coulomb phase, even if it is unobservable, has to be subtracted carefully in a numerical 
analysis in order to avoid that interesting finite terms are subtracted as well. Also, the 
branch point appears to be quite soft: a power {p — p*)'^/^ (p being a generic source 
parameter with p* its critical value) behaviour which has to be unraveled from leading 
analytic pieces (a constant and a linear term). 

Fortunately, a trick can be found in order to avoid the IR-singularity. Consider the 
partial derivative of the on-shell action wrt a particular parameter of the problem, e.g. 
the total cm energy. In general, this derivative will receive two contributions: one from 
the explicit dependence of the action from that parameter; and one from the implicit 
dependence of the solution (the fields) upon the parameter. However, this second contri- 
bution is also proportional to the variation of the action with respect to the fields which, 
by definition of the classical equations, vanishes on shell. We are thus led to: 

'-'P J \^P / fixed fields 

We can use this general strategy to our advantage in the case at hand by considering 
derivatives of ^ and, as parameter p, the overall cm energy B? while keeping the shapes 
of the two sources, Ri{t)/R fixed. Since: 



Gs 



J [ R^ R^p ^^ > I'y^^ 27rRpJ V 27cRp 



(30) 



we see that only the first term contributes to the abovementioned partial derivative so 

that: 

djA/Gs) 1 f 

^^^ = ^yrft(l-p) , (31) 

which is now perfectly IR-finite. 

Interestingly, the above derivative of the action is directly related to the average mul- 
tiplicity of emitted gravitons. Using the effective form of the complex action introduced 
in yj one readily obtains: 

which shows that the total multiplicity is in fact IR-finite in the axisymmetric case (recall 
that the IR-sensitive polarization is not produced in this case). 

We are now ready to analyze numerically the behaviour of the on shell action in the 
vicinity of the critical point. In Fig. [5] we show R'^dji2[A/Gs), as obtained from the 
integral fl3T|) in a range of i?'s. At first sight one might suspect that (N) is infinite exactly 
at i?c- A careful analysis shows that it is instead finite although it approaches the critical 
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point with an infinite derivative. The fit, shown by a sohd fine, is consistent with a square 
root branch point, 



(N) = Co + ci{Rc - R)-^ , 
in the vicinity of Re- This confirms the vahdity of the expansion 

A{R) =Ao + Ai{R, - Rf + A2{Rc - Rf , 



(33) 



(34) 



with 6=1 and c = 3/2. Such a behaviour, first found in [Ij and later confirmed in [7], 
could be possibly related [1] to a Choptuik-like [lOj dependence of the mass of the barely 
formed black hole from [R — Rc). In that case our result would be interpreted as meaning 
that Mbh ~ (-R — Rc)^, with 7 = 3/4, an exponent which is about twice the one found 
by Choptuik [TO] for the spherical collapse of some fiuids. It will be very interesting to 
see what values of 7 will come out of future numerical studies of the collapse of extended 
sources of the type discussed here. 



VI 




0.44 



0.445 



0.465 



0.47 



Figure 5: The total multiplicity of emitted gravitons (points) and the best fit: 0.138 — 
0.46(i?c — -R)°'^^'^. A fit with the fixed power 1/2 is marginally worse. 

One can also check directly the behaviour of the full, suitably regularized, action (1271) 
close to the transition point. We have used two regularizations: 1) the integral was cutoff 
at some large value tmax, defining the cutoff action Acut, and 2) we have introduced a 
counterterm, 1/ pit) — > ^/ pif) — 1/^) which subtracts the infrared divergence and results 
in the IR and UV-finite expression Agub- Figure [6] shows these actions together with fits 
like (IMl) and two a priori unconstrained powers. 

Comparing the two plots one sees immediately that the cutoff action is much more 
linear than the subtracted one. What happens is that A^ut is obviously much bigger (cf. 
the scales of both figures) and the nonlinear effect of the last term is relatively weaker. 
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0.36 0.38 0.4 0.42 0.44 0.46 



-0.02 - 



.36 0.38 0.4 0.42 0.44 0.46 



Figure 6: Fitting the critical behaviour of the cutoff (upper, tr, 
actions. 



50) and the subtracted 



On the other hand, in Asub the big "background" is already subtracted exposing better 
the last two terms. This observation is also confirmed by our fits. Simple Mathematica 
fitting routines were unable to disentangle two different powers in Acut, while fitting Ag^b 
gave b = 1.07 ± 0.12 and c = 1.34 ± 0.22 in good agreement with 6=1 and c = 1.5 
predicted in pj. As for Acut, we were glad to find that a slightly less general form. 



A{R) = Ao + Ai{Rc - Rf + A2{Rc - R) 



fe+l/2 



(35) 



works well, cf. Figl6] (upper part). Finally, we have also checked the cutoff dependence of 
these parameters, see Table [H 
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l^max 


b 


Ao 


Ai/Ao 


A2/A1 


20 


1.044 


0.826 


-6.149 


-0.894 


50 


1.045 


1.034 


-6.096 


-0.883 


100 


1.046 


1.189 


-6.077 


-0.877 


300 


1.049 


1.433 


-6.060 


-0.871 



Table 1: Parameters of the fit (!35l) for a range of cut-offs, tmax 

As expected the power b is cutoff independent and close to the predicted value 6 = 1. 
The amplitudes Ai show a mild, possibly logarithmic, dependence on tmax- Interestingly, 
the ratios of these amplitudes seem to depend very weakly on t^ax suggesting a possible 
factorization of the IR divergent contributions. 

6 Summary and future prospects 

In this paper we have continued the ambitious program initiated in ACV [Ij by focussing 
on the case of axisymmetric collisions. Both the ACV approach and the collapse criteria 
based on the explicit construction of minimal CTS [5], [6] greatly simplify in this case. In 
the former the axial symmetry removes the infrared problem and reduces the equations 
to ordinary differential equations, often allowing for quasi-analytical solutions. The CTS 
criteria also take a much simpler form, as exemplified by eq. flT4l) . 

In spite of its simplicity, this case becomes extremely rich if one generalizes the problem 
of particle-particle scattering to the one of the central collisions of two extended sources. 
The latter problem depends, in principle, on an infinite number of parameters (describing 
the energy profile of each beam) and thus allows to study, in principle, how the information 
about the initial state is transmitted to the final state. 

In this first paper we have concentrated our attention on the relevant system of equa- 
tions in position space, on the determination of the critical points/surfaces, and on the 
comparison with CTS-based results. Our results have confirmed a rather striking con- 
nection between CTS collapse criteria and the existence of critical points (in general 
hypersurfaces) in the parameter space describing the two incoming sources. Not only 
could we prove a general theorem that the CTS criterion of [6| implies the absence of 
real regular solutions of the ACV equations; we have also found that the bounds on the 
critical values based on [6] are typically only a factor two off the ACV predictions. 

We have also investigated the behaviour of the on-shell action, hence of the S-matrix, 
as one approaches the critical points, confirming the robustness of a definite fractional 
power-law behaviour for the phase of the S-matrix. Within some assumptions this might 
be interpreted as a critical "Choptuik exponent" [10] for the behaviour of the mass of the 
black hole as one approaches the critical point from the collapse side. 

Unfortunately, techniques for going beyond the critical points are still rudimentary, 
while the interpretation of the extra absorption of the elastic amplitude as being due to 
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the opening of new (black-hole formation?) channels remains very speculative. A full 
understanding of the role of inelasticity (i.e. of graviton emission) even in the subcritical 
regime is still necessary and might shed light of the nature of the "phase transition" 
between a "full dispersion" regime and one in which a fraction of the incoming energy 
collapses. For all these questions a momentum-space approach (see again [T|) -rather 
than the one in position space adopted here- may turn out to be more suitable. We plan 
to present it in a near-future investigation. 
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